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Abstract 

The Hamiltonian of the Rabi model is considered. It is shown that the ground 
state energy of the Rabi Hamiltonian is simple for all values of the coupling 
strength, which implies the ground state energy does not cross other energy. 



1 Introduction 

Cavity quantum electrodynamics has supplied us with stronger interaction than the 
standard quantum electrodynamics (QED) does |HR014 iRBHOlj . Experimental physi- 
cists usually demonstrate the interaction by coupling a two-level atom with a one-mode 
light (i.e., single-mode laser) in a mirror cavity (i.e., a mirror resonator). The region 
that the strong interaction in cavity QED belongs to is called the strong coupling 
regime. At the dawn of the 21st century, the solid-state analogue of the strong interac- 
tion in a superconducting system was theoretically proposed in |MSS01t iMBOlj . and it 
has been experimentally demonstrated in |CBSNHM04|IF(;BBLBW08llWSBFHMK(;S04| . 
That is, the atom, the light, and the mirror resonator in cavity QED are respectively 
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replaced by an artificial atom, a microwave, and a microwave resonator on a supercon- 
ducting circuit. Here, the artificial atom is made by using a superconducting circuit 
based on some Josephson junctions. This replaced cavity QED is the so-called circuit 
QED. The circuit QED has been intensifying the coupling strength so that its region 
is beyond the strong coupling regime. This amazing region of the coupling strength 
between the artificial atom and the light is called the ultra-strong coupling regime in cir- 
cuit QED |DGS07l IGXHSBSLCTLH091 IFLMGSHMTOl INDHMHSGZHSMGTO] . Then, 
experimental physicists have found some differences in physical phenomena between the 
two coupling regimes jFLMGSHMlOl INDHMHSGZHSMG10] . As one of the striking 
differences, there is the following. In the strong coupling regime as well as in the weak 
coupling regime, the Jaynes-Cummings (JC) model is useful to explain the experimen- 
tal results jHROll IFGBBLBWl)8] . The Hamiltonian of the JC model is obtained by 
applying the so-called rotating wave approximation (RWA) to the Rabi Hamiltonian. 
On the other hand, in the ultra-strong coupling regime, the JC model does not work, 
and thus, we need a help of the Rabi model |FLMGSHMlUl INDHMHSGZHSMGlOj . 
The current cutting-edge technology of circuit QED is begining to show us the divi- 
sion between the two coupling regimes concretely. We are interested in how physics 
determines this division. 



To see a difference between the Rabi model and the JC model, in this letter we 
pay particular attention to the energy level crossing. Recently, Braak |Brallj had 
given a mathematically intriguing expressions of the eigenenergies of the Rabi model. 
Then, the following questions arise and are problems of interest to us: (i) are there any 
energy level crossings among them? If so, (ii) how do they take place? As is shown in 
[Hir09at lHir09b| . the ground state energy of the JC model consists of many energy level 
crossings as the coupling strength grows larger and larger. Namely, for the JC model 
the quantum phase transition in Rey's sense |Rey09| takes place (see FigH] below). For 
the details on energy level crossing and quantum phase transition, see |Sach 99j. On 
the other hand, as in our numerical computation in Fig|2j we can conjecture that the 
ground state energy of the Rabi model has no energy level crossing. In this letter 
we prove this fact with the functional-integral method [HLQZl IHIL12| [HHL12J as a 
corollary of the fact stating that the ground state energy of the Rabi model is simple 
(i.e., the ground state is unique). It reveals us that it is in the ultra-strong coupling 
regime of circuit QED that there is a big qualitative difference as well as quantitative 
one between the Rabi model and the JC model. 
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2 Rabi model 
2.1 Definition 

Let (T x ,<7 y , a z be the 2x2 Pauli matrices: 
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In this letter we adopt the natural unit: H — 1. The renormarized Hamiltonian of the 
Rabi model is defined as a self-adjoint operator by 

#Rabi = cr z A + uja ] a + ga x (a + a 1 ) (2.2) 

on C 2 <S) L 2 (r). Here A > and u > are respectively the atom transition frequency 
and the cavity resonance frequency, g ER stands for a coupling constant, and a and a) 
denote the single mode bose annihilation and creation operators satisfying [a, a*] = 1. 
It is given by a = + y/ux)/ 'y/2 and a) = (— + y/ux)/y/2. We are interested 

in studying spectral properties of eigenenergies of i^Rabi, in particular crossing of the 
ground state energy. 

The absence of crossing can be derived from the simplicity of the ground state 
energy of i?R a bi- We will construct a path integral representation of e _< ^ Rabi to show 
that the ground state energy is simple. This is a minor modification of recent papers 
[HLOT] IHIL12[ IHHL12] , where the Feynman-Kac type formula with spin is established. 
In particular the spin-boson model is studied by path measure in |HHL12] and we can 
apply it in this paper since the Rabi model can be regarded as the single mode photon 
version of the spin-boson model. 



2.2 Two conjectures 

Let us here consider the Rabi Hamilonian before the renormalization: 

H R = H Rahi + 1 = a z A + u (^a ] a + + ga x (a + a ] ) . 

In this letter we follow the clasification proposed in [CRLGS10J, and define the ultra- 
strong coupling regime by the region in which the dimensionless coupling strength 
g/u > 0.1. 

Applying the RWA to H R , we have the JC Hasmiltonian: 

Hjc = ff 2 A + w (a^a + i J + g (a^a) + a + a) , 
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where a± = (a x ± ia y )/2. We denote by Ejq the ground state energy of Hjc- The 
JC model is a completely solvable model, and the eigenstate fi c (g) of H JC and its 
corresponding eigenenergy E'l c (g) are given for each v G Z in the following procedure: 
Let tp g (x) = (u/ixy^e^ 2 / 2 be the normalized eigenvector associated with the lowest 
eigenvalue, 0, of the harmonic oscillator ua^a = +oj 2 x 2 — oj). Then Fock states 

are defined by \n) = ^=(n™ a ^)v 9 g f° r the single mode photon, n = 0, 1,2, • • • , with 



|0) = ip g . For the spin ground state |— ) = 

Acr 2 , we define states | — ,n) and |+, n) by 
respectively Then, 



and the spin excited state |+) = 
n) = |— ) ® \n) and |+, n) = |+) 



of 

n), 



'^ c (s) = |-,o>, 

<f^ v \(j9) = cos(0 H (#))|+, \u\ - 1) + sin(0 H (#))|- u^O, 
V%\{9) = ~ sm(e H (g))\+, \u\ - 1) + cos(0 M ((/))|- i/ ^ 0, 

where H (#) = | tan" 1 ^^§J if 2A ^ w; H (#) = tt/4 if 2A = w, and 



g^i(g)=^H±V (2A 4 a ' )2 +g 2 H. 

According to |Hir09at IHir09b| , the remarkable finding for Ejq is the energy level cross- 
ings in the ultra-strong coupling regime: For each n — 0, 1, 2, • • • , there exists g n+ i > 
such that E^(g) and E^9 n+1 Jg) cross each other at g — g n +i, and 

£jc = E J _ c n (g) ifg< g n+1 , 

E JC = Ei c n (g) = E™(g) if g = g n+1 , 



-(n+l)» 

.£jc = £^ n+1) (#) if#># n+ i, 



provided 2A > u. See FigjTj In other words, as g gets large, there exists v g G Z_ 
such that i^jc = E^(g), and moreover, 2/ s is strictly decreasing and v g — >■ — oo as 
<7 — >■ oo. Namely, these energy level crossings take place and make the ground state 
energy Ejc as the envelop of El°(g), v = 0, —1, —2, • • • , in FigjlJ We also note that 
the ground-state entanglement [PZDS10J for the JC model. Namely, for instance, the 
ground state is a separable state for g < g 1; but it becomes an entangled state for 
g > 9i- The details on g n and u g are in |Hir09at IHir09b] . 

In Fig|2] there is a numerical computation of the energy levels of H^. It says that 
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Fi glire 1: Energy level crossings among E^ c (g), v = 0, —1, —2, ■ ■ ■ , of the JC Hamiltonian. Each color indicates 
individual index v of the energy E^p(g). (a) 2A = ui; (b) 2A = 3ui. 
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0.5 1 1.5 2 2.5 
strength g/co of interaction 

Fi gUre 2: Energy levels of for 2A = ui. Each color indicates the nth level of the energy for n = 0, 1,2, • • • from 
the bottom, where the 0th level energy means the ground state energy. 

(CI) there is no energy level crossing between the ground state energy and the 1st 
excited state energy; 

(C2) we may say that there are n energy level crossings between the 2n-th excited 
state energy and the (2n + l)-th excited state energy, n = 1, 2, • • • . 

In this letter, we will prove (CI). 

3 Results and proofs 

Before going to show the Feynman-Kac formula of e~ Habi , we prepare a probabilistic 
description of H^ahi- 

Let <j = (<j x , <j y , <j z ) be elements of SU(2). The rotation group in E 3 has an adjoint 
representation on SU(2). Let n G R 3 be a unit vector and 9 G [0, 27r). Thus we have 
e (i/2)8n-* satis fi es that e (i/2)en-a a ^ e -(i/2)e n -o = (ifo)^ w here R denotes 3x3 matrix 

representing the rotation around n with angle 9. In particular for n = (0, 1, 0) and 

9 = 7T/2, We have e m)en-a axe -(i/2)9n-a = Gz and e (i/2)9n-a aze -(i/2)en-a = 
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U = e (i7r/4)CT «. Then 



UH Rahi U 1 = uja ] a + ga z (a + a 1 ) - a x A. 



(3.1) 



Since <p g is strictly positive, we can define the unitary operator U g : L 2 (R) — >■ L 2 (]R, y>gd:r) 
by U g f = ^Pg 1 /- We set the probability measure (p 2 dx on K by d/i. Thus UH^uU" 1 
is transformed to the operator: 



UgUH^uU- 1 !/- 1 = - ( -— + ui- J + ga z ^x - a.A 



(3.2) 



in C 2 <S> L 2 (]R, d/i). Let us introduce Z 2 = {— 1,+1} to redefine the Hamiltonian (3.2) 
on a set of scalar functions. We identify C 2 ® L 2 (R, dfi) with 



= L 2 (R x z 2 , d/i) = < / = f(x, a) 



(3.3) 



[/-(*). 

transformed to the operator if: 



by C 2 <g) L 2 (K, d//) 9 



1— )■ f(x,a) G Jf. Thus under this identification (3.2) is 



1 



d 2 d\ 1 

+ ux— I + gV2uax > f(x, a) - Af(x, -a), a G Z 2 (3.4) 



2 y da; 2 ' da; 
in Jff. Thus we have the lemma below: 

Lemma 3.1 The operator iiRabi in C 2 £§> L 2 (R) is unitarily equivalent to H in 3t? 

In what follows we deal with H instead of H Rsihi . Let 

1 



h 



d 2 d 
dx 2 dx 



and (X t ) t >o be the Ornstein-Uhrenbeck process on some probability space (C, S,P X ) 
We have P X (X = x) = 1 



dfi(x)E P x [X t ] = 0, 
dfi(x)E P x [X t X s ] = 



-\t-s\uj 

2oo 



Here Eq [■ ■ ■] denotes the expectation with respect to a probability measure Q. The 
generator of X t is given by — h and 



(f,e~ th gU= ldv{x)E P * f(X )g(X, 
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The distribution p t (x,y) of X t under P x is given by 

Pt(x, y) = (fgix^Ktix, y)<p s (y) 

where K t (x,y) denotes the Mehler kernel given by 

1 



(3.5) 



K t (x,y) 



-2t\ 



( Axye' 1 - (x 2 + y 2 )(l + e~ 2t ) 

exp l 



-2t\ 



See e.g., |LHBlll 3.10.4] for the detail of Ornstein-Uhrenbeck processes and harmonic 
oscillators. In order to show the spin part by a path measure we introduce a Poisson 
process. Let (N t ) t > be a Poisson process on some probability space (C, £', v) with 
unit intensity, i.e., 



t n 

E u [l Nt =n] = 

n! 



We define a t = (-l) Nt a, a G Z 2 , for t > 0. Let 

^ d/ji(x)EpxE v [■ • 
o-ez 2 



n > 0. 



E 



(3.6) 



Theorem 3.2 It follows that 

(A>0) (/,e-^)^ = e*E 



(A = 0) (/,e-*^) 



f(X , a )g(X t} at ) e -9V^J a s x ad s A N t 
t Nt=0 f(X ,a)g(X t , a ) e -9^S*x 3ds 



(3.7) 
(3.8) 



Proof: Let A > 0. By a minor modification of [HIL12, Theorem 5.10] we can see that 



(f,e~ tH g)j? = e*E [f(X ,a )g(X t ,a t ) 



gy/2L0 f a a X a ds „/ log A<£/V a 



(3.9) 



Here / f(N s )dN t = f( N r)- Since e^ ogAd7Vs = e logAiVt = A Nt , (3.7) follows. 

"^0 <r> AT . -A AT 



r,N r+ ^N, 

In the case of A = only the set N t = contributes to the path integral. 



Then 



(/, e~ tH gU = e'E f(X , a )g(X t , a t ) e -*>^ 



N t =0 



Then (3.8) follows. 



(3.10) 
qed 



Let E = info- (if). 



Corollary 3.3 It follows that dimker(H — E) = 1, i.e., the ground state of i^Rabi is 
unique. 
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Proof: Let /, g > but not identically zero. Then for sufficiently small e > 0, we see 
that both £lf = {(x,a) ERx Z 2 \f(x,a) > e} and £l g = {(x,a) ERx Z 2 \g(x,a) > e} 



have positive measures. We have by (3.7) 



Since Qf is a subset of MX Z 2 , we have Qf = Uo-ez 2 (^/' ° r )- Thus either fijr or f2j (c M) 
have at least a positive measure. Suppose that has a positive measure. Similarly we 
see that fl g = Uo-ez 2 (^9' a ) an< ^ su PP ose that f2+ is a positive measure. Let Q be the 
set of paths starting from the inside of , +) and arriving at the inside of +)• 
We see that 



E[Io]=E I Qi (X )l a +(X t )l 



Nt =even 



By using the distribution p t of X t we have 



E 



00 f2n r r 



dx dyip g (x)K t (x,y)ip g (y) > 0. 



Hence we conclude that has a positive measure and 



Cf,e-* H <7)>ee*E 



> 0. 



Thus e tH is a positivity improving operator. Thus dimker(i/ — £') = ! follows from 



the Perron-Frobenius theorem. 



qed 



Corollary 3.4 The ground state energy of i?R a bi has no crossing for all the values of 
g and A. 



4 Conclusion 

In this letter we have proved the first conjecture (CI) that the numerical computation 
predicts, while the JC model has many energy level crossings for the ground state 
energy in the ultra-strong coupling regime of circuit QED though it has no energy 
level crossing in the weak and strong coupling regimes. It shows that it is in the ultra- 
strong coupling regime that there is a big qualitative difference as well as quantitative 
one between the JC model and the Rabi model. 
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